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Abstract 

We give a brief overview of nuclear parton distributions. First, the EMC effect is dis- 
cussed together with possible interpretations such as nuclear binding and Q 2 rescaling. 
Next, we explain shadowing descriptions: vector- meson-dominance-type and parton- 
recombination models. Nuclear dependence of Q 2 evolution should be interesting in 
testing whether or not DGLAP equations could be applied to nuclear structure func- 
tions. Status of nuclear sea-quark and gluon distributions is discussed. The structure 
function b±, which will be measured at HERMES and possibly at ELFE, could shed 
light on a new aspect of high-energy spin physics. 

1 Introduction 

From the early SLAC experiments on electron scattering, nuclear targets have been 
used just as a tool for measuring nucleon structure functions. Nuclear corrections are 
not taken seriously in 1970's except for apparent nucleon-Fermi-motion corrections. 
Nucleons are bound in a nucleus with average binding energy 8 MeV and typical 
momentum 200 MeV. These are fairly small compared with reaction energy of several 
GeV or more in electron or muon scattering. For this reason, nuclear environment was 
not expected to alter nucleon properties significantly. Even though signature of nuclear 
modification in the F2 structure function is found in the SLAC data in 1970's, they are 
not studied extensively. On the other hand, the European Muon Collaboration (EMC) 
took the topic rather seriously and published the first paper in 1983 ||. The nuclear 
modification of F2 is now known as the (old) EMC effect. 

In the mid 1980's, there are many theoretical publications on this topic, in particular 
on medium and large x physics. We cannot list all of the ideas in this paper. A 
conservative one is to interpret it in term of nuclear binding and Fermi motion || . On 
the other hand, there is another extreme model, Q 2 rescaling [0]. In contrast to the 
binding model in which bound nucleon F 2 is the same as the free one, it is modified 
due to nuclear environment. The details of these models are discussed in section |2|. 

From the late 1980's, nuclear F2 structure functions in the shadowing region (x < 
0.1) are investigated in detail experimentally and theoretically. In particular, accurate 
New Muon Collaboration (NMC) data make it possible to test theoretical shadowing 
models. The traditional idea of describing the shadowing is a vector-meson-dominance 
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ing has been used for probing internal struc- 
ture of nucleons or nuclei. Its cross section 
is written in terms of two structure functions 
Fi and F 2 . F\ is related to the transverse 
cross section for the virtual photon, and F 2 
is to the transverse and longitudinal ones. 

If higher-order a s corrections are neglected, 
Figure 2.1: Nuclear modification F? a /Ff. , , , , ^ , ,. 

° - they are related by the uallan-Gross relation 

2xF\ = F 2 . Therefore, we discuss only the F 2 structure function. In parton model, 

the cross section is described by the incoherent summation of scattering cross sections 

from individual quarks: 

F 2 (x) = X)e?x[gi(x)+ft(x)] . (2.1) 

i 

We discuss how the structure function F 2 is modified in nuclei. The details of the 
modification is studied experimentally for several nuclei. For example, SLAC, EMC, 



and NMC data || for the calcium-deuteron F 2 ratio are shown in Fig. |2.1] . In the 
medium x region [x ~0.5), the ratios are approximately 10% smaller than unity, and it 
tends to increase at large x(> 0.8). In the region x ~ 0.2, the modification is positive, 
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and the phenomenon is called "antishadowing" in contrast with shadowing at small x. 
The shadowing part is discussed in section 

In describing nuclear structure functions, convolution formalism is usually used. 
A parton distribution in a nucleus is given by a hadron (nucleon, pion, and etc.) 
momentum distribution convoluted with a parton distribution in the hadron. In terms 
of light-cone momenta, it is written as 

fa/A(x A ,Q 2 ) =Y1 I d y A f a/T (x A /y A ,Q 2 )f T/A (y A ) , (2.2) 

where x A is defined by x A = Q 2 /2M A v and y A is the light-cone momentum fraction 
Va = Pt/Pa- m order to explain the modification of the quark distribution in a nucleus 
[fa/A(x A , Q 2 )]i we ma y ascribe it to modification of the quark distribution in a hadron 
[fa/T(x A /y A )} or to modification of the hadron distribution in the nucleus [/t/aCz/a)]- 

A standard approach in describing nuclei in the low-energy region is to start from 
information on nucleon-nucleon interactions, which are obtained by NN scattering ex- 
periments and deuteron properties. Then, the knowledge is applied to nuclear structure 
by taking into account many-body corrections. In this sense, it is a "proper" approach 
to start from an assumption that the structure function F 2 of a bound nucleon is the 
same as the free nucleon F 2 . This way of interpreting the EMC effect is a nuclear bind- 
ing model [0]. A nuclear tensor, which appears in calculating a lepton-nucleus cross 

section, „ gi ven by *. — „ WW - The S pe Ctral 

function S(p) is the nucleon momentum distribution in the nucleus. If we consider the 
simplest case of a shell model, it is given by single particle wave functions and binding 
energies of nucleons: S(p) = EI ( Mp)| 2( Kpo — — e»). Recoil energy is neglected in 

i 

the expression. From these equations, the nuclear F 2 is 

F 2 A (x, Q 2 ) = E / dzfiWFfix/z, Q 2 ) , (2.3) 

i 

where x is the Bjorken variable x = Q 2 /2m^i>, and fi{z) is the light-cone momentum 
distribution of the nucleon i: 

f t (z) = [d 3 pz5(z- P^L) MP)\ 2 ■ (2.4) 

The momentum fraction z is given by z = p ■ q/m^v = 1 — |ei|/mjv + P • q/^n^ ~ 
1.00 — 0.02 ± 0.2 for a typical nucleon. So, the function f(z) is peaked at 0.98. If 
we replace f(z) by a delta function 5(z — 0.98) for simplicity, the nuclear F 2 becomes 
F 2 A (x,Q 2 ) w F 2 N (x/0.98,Q 2 ). For example at x=0.60, we have Ff{x = 0.61)/F 2 JV (x = 
0.60) = 0.88 and the 10% EMC effect can be explained at medium x. 

On the other hand, there exist other extreme interpretations in terms of modifica- 
tion of nucleon itself. A well known idea is the Q 2 rescaling model. It seems that the 
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nucleoli modification makes sense if we consider that the average nucleon separation in 
a nucleus is about 2 fm, which is almost equal to the nucleon diameter. Therefore, it is 
no wonder that multiquark systems other than the nucleon are formed. If such multi- 
quark hadrons are created, a confinement size for quarks should be changed. Then the 
quark momentum distributions are modified according to the size change. Using this 
kind of simple picture at a hadronic scale /i 2 , we calculate distributions at large Q 2 , 
where experimental data are taken, by evolving distributions from fi 2 . We find that a 
nuclear F 2 is related to the nucleon one simply by scaling Q 2 \ 

F 2 A (x,Q 2 )=F 2 N (x,UQ 2 ) , (2.5) 

with the rescaling parameter £4 = (\\/ 'A 2 v ) Qs ^a)/ Qs (Q 2 ) ? where A is the confinement 
radius for a quark. 

The binding and rescaling models seem very different interpretations. How shall 
we understand the situation? A possible explanation is in terms of factorization 
scale independence Q. For example in the operator- product-expansion case, the 
renormalization scale /1 separates long-distance physics and short-distance one. How- 
ever, final results should not depend on the arbitrary human factor \i. Because 
the bound structure function and the nucleon momentum distribution are not sep- 
arately observables, the same discussion could be valid in our case. The nuclear 
structure function, which is a physics observable, should not depend on a factor- 
ization scale n which separates it into the quark-distribution part and the hadron- 
distribution one. Choosing a nuclear-dependent scale /1 = /i^, we obtain the rescaling 
model. On the other hand, the binding model corresponds to the nuclear-independent 
scale. Using this factorization-scale independence, we may relate the two different de- 
scriptions. Moments of a nonsinglet distribution are written as M^J A = M% N M^' A 
by Mellin-transformation properties. As a simple example of the binding model, a 
Fermi-gas model is employed. Relations between the two descriptions are found as 
e/m N = -^ s k a /2(3 and k 2 F /m 2 N = 5(2 7 f 5 - ^ s )k a /2(3 + 0{f/m 2 N ), where k a is 
defined by 1 — k a = ® s (£,aQ 2 ) / &s{Q 2 ) ■ The quantities e and kp are the average bind- 
ing energy and the Fermi momentum, and 7^ is the nonsinglet anomalous dimension. 
There exists a correspondence between the binding model and the rescaling, so that 
we may view the rescaling as an effective description of nuclear medium effects. In 
fact, it is shown in a simple quark model that effects of nucleon interactions could be 
effectively described by a nucleon size change in the nuclear medium || , although it is 
perhaps too simple to explain various nuclear phenomena. Comparison of the rescaling 
results with the data is discussed in section 19. 



3 Nuclear shadowing 

Nuclear modification of F2 at small x is negative as shown in Fig. 1271], and it is known 



as shadowing phenomena. A traditional way of describing the shadowing is in terms 
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of vector-meson-dominance (VMD) model. The virtual photon transforms into vector 
meson states, which interact with the target. Because a propagation length of the 
vector mesons is given by A = l/\Ey — E^\ «0.2/x fm, it exceeds 2 fm at x <0.1 and 
multiple scattering occurs. The vector mesons interact predominantly in the surface 
region, and internal nucleons are "shadowed" by the surfaces ones. This description is 
not enough for explaining observed small Q 2 dependence. So various extensions of this 
model are studied. For example, including qq continuum, we have the nuclear F 2 JF| 

-g/^ ^^U , , (3.1) 
2 V ,yi ' ir J (M 2 + Q 2 ) 2 V 7 

where ay a is the vector- meson nucleus interaction cross section, and II (M 2 ) is the spec- 
tral function II (M 2 ) = cr(e + e~ — ► hadrons)/a(e + e~ —>■ fi + fi~) = vector mesons + 
qq continuum. 

Another interpretation of shadowing is in term of Pomeron exchange ||. In the 
case of diffractive scattering, the target is thought to be remain intact so that only vac- 
uum quantum number, "Pomeron" , could be exchanged. The Pomeron (V) structure 
function is defined by the diffractive cross section: F 2 j> = <3 2 /(47r 2 a)<T 7 *p. Pomeron 
contribution to the nuclear F 2 from double diffractive scattering is given by a con- 
volution of the Pomeron F 2 with its light-cone momentum distribution: 5F 2 (x) = 

The — , i3 the _ m fra e t ,oe eaeried h Y t he P 0merOn 

U = k ■ q/p ■ q. The VMD contribution is compared with the Pomeron result in |J. 
Both shadowing results are of the same order of magnitude at small x (~0.01) and at 
Q 2 =4 GeV 2 . 

The above shadowing ideas are in the target rest frame. There is an explanation 
in an infinite momentum frame in terms of parton recombinations (|Tl], [TT|] . Both 
descriptions are supposedly equivalent; however, there is no explicit proof at this stage. 
The recombinations are parton interactions from different nucleons. They become 
dominant at small x with the following reason. The average longitudinal nucleon 
separation in a Lorentz contracted nucleus is L = (2 /m)M^/ 'Pa = (2 fm)m N /p N . 
The longitudinal localization size of a parton with momentum xpj^ is AL = l/xpjy. 
In the x region x < 0.1, AL becomes larger than L, so that partons from different 
nucleons could interact significantly. For example, &pip 2 —>■ Pz{x) recombination effect 
on the P3 parton distribution is 

9A 1 / 3 a s r 

Ap 3 (x) = s / dx l dx2Pi(xi)p 2 (x2)T PlP2 ^ P3 (x 1 ,X2,x)6(x -x x - x 2 ) , (3.2) 

where the factor A 1 / 3 is the number of nucleons in the longitudinal direction, p±p 2 — ► p% 
cross section is proportional to a s /Q 2 , and Rq is given by Ro ~1 fm. The parton-fusion 
function r pi p 2 ^p 3 indicates a probability of the fusion process pip 2 — > p%. Because it 
is opposite to the splitting, r piP2 ^p 3 is related to the Altarelli-Parisi splitting func- 
tion by T pip2 ^ P3 (x 1 ,x 2 ,x 3 ) = (x 1 x 2 /xl)P pi ^ m (x 1 /x 3 )C PlP2 ^ P3 , where C P1P2 ^ P3 is the 
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show the nuclear modification F 2 jF 2 in a 
model with the recombination and the rescal- 
ing mechanisms |TTJ. Both contributions are 
calculated at Qq, then obtained distributions 
are evolved to the ones at larger Q 2 . The ini- 
tial point Ql is considered parameter in 
the model, and it is fixed by fitting the NMC 
shadowing data for the calcium (Ql=0.8 GeV 2 ) 
The parton distributions at Q 2 =5 GeV 2 are 
calculated by the DGLAP evolution equations, 
and they are compared with the data in Fig. |3.1| . Although Ql is an adjustable param- 
eter, the model can explain the nuclear modification from very small x to large x. It is 
interesting to note that the antishadowing part at x ^0.2 is explained by competition 
between the recombination and the rescaling effects and that the "Fermi motion" part 
is described by the quark-gluon recombinations. 



Figure 3.1: Comparison of theoretical re- 
sults with the data F§ a jF§ ' . 
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Table 3.1: Status of nuclear parton distributions. 

Current status of nuclear parton distributions are summarized in Table |3.1| . The 
F 2 structure function at medium (small) x is essentially the valence (sea) quark dis- 
tribution. So we have good information on the valence-quark distribution at medium 
x and the sea-quark one at small x in several nuclei. The sea-quark distribution is 



also investigated in Drell-Yan processes. The Fermilab Drell-Yan experiment [12] did 



not measure the shadowing region as shown in Fig. 3.2. We hope that future mea- 



surements at Fermilab and at RHIC probe the shadowing region at very small x. The 



iron data in Fig. |372| are often quoted for concluding that there is little modification 
in the sea. However, it is not very obvious by looking at other nuclear data. Accurate 
A dependent data, as well as theoretical studies, are necessary for finding the details 
of sea-quark behavior. Because the sea-quark distribution dominates F 2 at small x, 
behavior of valence quarks at small x is not known. However, it could be an inter- 



esting topic in testing various shadowing models |13| . The valence-quark shadowing 
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Figure 3.3: Gluon distribution in muon- 
induced 3/ijj production. 



4 Nuclear dependence of Q 2 evolution 



It is known that the Bjorken scaling, which indicates structure functions are indepen- 
dent of Q 2 , works approximately. Therefore, most studies on the nuclear modification 
are focused on the x dependence. It is difficult to find small nuclear effects on the scaling 
violation. However, it became possible to measure such small effects recently The 
NMC obtains Q 2 evolution differences in tin and carbon nuclei: diF^/F^ )/d(\n.Q 2 ). 
It is an interesting topic to investigate whether or not the DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) equations could be applied to nuclear structure functions. In 
particular, the parton recombinations are predicted to produce higher-twist effects in 
the evolution 



WX 15 



= £j X ''/,,. (j) Qjfat) + P„\jjJ g(y,t) 

( , a s A 1/3 \ 

+ recombination terms oc — — — , (4.1J 

V Q J 

( a s A 1/3 \ 

+ recombination terms oc s — , (4.2) 

V Q J 



where the variable t is defined by t = — (2//5 ) \n[a s (Q 2 ) /a s (Ql)}. The first two terms in 
Eqs. (4.1) and (4.2) describe the process that a parton pj with the nucleon's momentum 
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nuclear dependence in Eqs. (4.1) and (4.2). 
First, parton x-distributions are modified in 
a nucleus so that input distributions on the 
right-hand sides are nuclear dependent. Sec- 
ond, the additional recombination terms have 
A 1 / 3 nuclear dependence. We calculate these 
effects on the nuclear difference in the Q 2 evo- 
Figure 4.1: Difference of Q 2 evolution in tin lution: d{F$ n / F°) / d{\nQ 2 ) [TjJ. As initial 
and carbon nuclei: d(F$ n / F?)/d(hiQ 2 ). distributions at Q%, we use those in the re- 
combination model with the rescaling in sec- 
tion |3|. They are evolved to Q 2 =5 GeV 2 by leading-order (LO) DGLAP, next-to- 
leading-order (NLO) DGLAP, and parton recombination (PR) equations with the help 



of the computer program bfl.fort77 in Ref. fj"5(|. The results are compared with the 



NMC data [|J in Fig. |4.1| . The DGLAP evolution curves agree roughly with the experi- 
mental data, so that the major nuclear dependence comes from the nuclear modification 
in the x distributions. The PR evolution curve disagrees with the data; however, it 
does not mean that the recombination mechanism is ruled out. There are a few uncer- 
tain factors. For example, the input higher- dimensional gluon distribution is essentially 
not known. Further refinement of the recombination contributions to the evolution is 
necessary. It is nonetheless interesting in the sense that higher-twist effects could be 
found in studying the nuclear Q 2 evolution. 



5 Spin structure of spin-one hadrons 

Spin structure of spin-1/2 hadrons has been investigated fairly well. In addition, there 
is a good possibility to measure a new spin structure function for spin-one hadrons in 
the near future. As a realistic spin-one target, we have the deuteron. Studying the 
deuteron spin structure at high energies, we gain an insight into a new aspect of spin 
physics fll7j| . There are four new structure functions for spin-one hadrons: bi, b 2 , b 3 , and 
6 4 . In the Bjorken scaling limit, the only relevant structure function is bi or equivalently 
62/2X. For measuring 6 1; the electron or muon does not have to be polarized but we need 
polarized deuteron. The b\ structure function is proportional to the tensor combination 
of the polarized cross sections: b\ oc da(0) — [da(+l) + da(— 1)]/2. In parton model, it 
is given by 
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where Sqi = [q®— (g i +1 +g i " 1 )/2]/2 = \\qi(x) — qf l (x)] with the flavor- i quark distribution 
q^(x) in z component of the target spin (H). Because quark spin does not appear in 
the above expression, b\ probes very different spin structure. 

We discuss a sum rule for b% ||18|| . The integral of Eq. ( |5.1| ) over x is related to 
a tensor combination of elastic amplitudes by using the parton picture in an infinite 
momentum frame. Next, the elastic helicity amplitudes are written by charge and 
quadrupole form factors of the spin-one hadron. In the case of tensor combination, the 
charge form factor cancels out and the remaining term is the quadrupole one: 



J dx bi(x) 



lim 

t->o 



(5.2) 



where the sea-quark tensor polarization is defined by 5Q sea = J dx[85u + 25d + 5s + 5s], 
and Fqit) is the quadrupole form factor in the unit of e/M 2 . Equation (|5.2|) is analogous 
to the Gottfried integral J dx \F%(x) - F^(x)} = 1/3 + (2/3) J dx[u{x) - d(x)}. As the 
Gottfried sum 1/3 is obtained by assuming the flavor symmetry u = d, the b\ sum 

5 t 



( dx bAx) = lim 77^F Q (t) = 

J v ; t->o 3 4M 2 QK ' 



(5.3) 



is obtained if there is no sea-quark tensor polarization. Because of this assumption on 
the tensor polarization, the b\ sum rule is considered as one of the sum rules in a naive 
parton model. It is similar to the Gottfried sum rule in this sense. At this stage, all 

model calculations for bi{x) satisfy the sum rule J dx b\(x) = 0. As the breaking of the 

Gottfried sum rule became an interesting topic recently, it is worth while investigating 
a possible mechanism to produce the tensor polarization 5Q sea , which breaks the sum 
rule. 

Even though the sum- rule value is expected to be zero for the b\, it does not mean 
that bi itself is zero. In fact, bi can be negative in a certain x region. For example, the 
deuteron b\ is given by a helicity amplitude for the nucleon convoluted with a light-like 
momentum distribution of the nucleon ]17| |19|| : 



f ^F x (x/y) I d 
Jx y J 

x (3 cos 2 9 - 1) 



V 



sin a cos au s (p)ud(p) + ——- sin 2 au 2 d {p) 



pcos9 

1 + — + 



P 



AM 2 
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y - 



16tt 

p cos 9 + E{p) 
M 



(5.4) 



where sina is the D-state admixture, and u Sj d are S and D-wave deuteron wave func- 
tions. The first term is due to the S-D interference and the second to the D-state. This 
is just an example of nuclear structure aspects. Measured deuteron b% does not have 
to agree with the above estimate. If a deviation from nuclear contributions is found, 
b\ provides important clues to physics of non-nucleonic components in nuclei and to 
new tensor structure on the parton level. 
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